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COOLING OF GAS TURBINES 
VIII - THEORETICAL TEMPERATURE DISTRIBUTIONS 
THROUGH GAS TURBINE WITH SPECIAL BLADES 
AND COOLING FINS ON THE RIM 
By W. Byron Brown, and John N. B. Livingood 


SUMMARY 

A theoretical analysis of the radial -temperature distribu- 
tion through a turbine with specially designed turbine blades was 
made for a turbine with and without cooling fins on the rim. The 
effect on blade life and effective gas temperature of the addi- 
tion of cooling fins on the rim was investigated for selected 
turbine operating conditions. The accuracy' of the radial- 
temperature distribution was then determined by finding a two- 
dimensional temperature distribution through the turbine rim and 
rotor for the turbine with and without cooling fins on the rim. 

The two-dimensional distribution was obtained by consideration 
of radial- and axial -temperature gradients and by application of 
the relaxation method to the resulting differential equation. A 
three-dimensional temperature distribution for a section of the 
turbine rim near a blade root was then determined by use of three- 
dimensional relaxation. 

The results showed that, for the selected conditions, the 
addition of cooling fins on the rim permitted only a slight Increase 
in the effective gas temperature. The two-dimensional investigation 
proved the radial distribution to be sufficiently accurate for most 
applications. The three-dimensional study indicated the type of 
temperature gradients to be expected in a section of the rim sur- 
rounding a blade root and may be applied to determinations of rim 
thermal stresses. 


INTRODUCTION 

High -temperature materials now available for gas turbines 
limit turbine -inlet temperatures to about 1500° F; consequently, 
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some method of turbine cooling is necessary if inlet temperatures 
above this value are to be used. Investigations have been con- 
ducted at the NACA Cleveland laboratory (references 1 to S) to 
determine the relative merits of indirectly cooling turbine blades 
(by conducting heat away from the blades) and of directly cooling 
them (by passing air or liquids through hollow passages in the 
blades). Increases in gas temperatures of the order of 100° to 
200° -F were found possible for indirect cooling; considerably 
greater increases in gas temperatures were made possible by direct 
cooling with either air or water as the coolant. 

In order to apply direct cooling using either air or water as 
the coolant, the turbine blades must be designed with hollow pas- 
sages. The manufacturing problems involved in the design of hollow 
blades, however, appear to limit blade forms as compared with solid 
blades. Seme turbine blades for indirect cooling have been designed 
and tested; one of these is a blade with a long, thin trailing sec- 
tion and with a high degree of twist. 

The results of a theoretical investigation of indirect cooling 
for a turbine with this special blade are presented herein. Eadial- 
temperature distributions were obtained before and after the addi- 
tion of cooling fins on the rim and the increases in blade life and 
effective gas temperature resulting from the addition of rim fins 
were calculated for a chosen set of turbine operating conditions. 

The radial- temperature distribution was then used as a trial solu- 
tion and a two-dimensional temperature distribution through the 
turbine rim . and rotor was obtained by considering additional tern-, 
perature gradients in the axial direction. A three-dimensional 
relaxation was then carried out in a section of the rim surround - 
;Lng a blade root. 


SYMBOLS 


A 

B 

C 

D 

E 

F 

G 


Tie following symbols are used in this report : 



average cross-sectional area, 


(ft) 


(°F) 

> integration (°F) 
| constants 


(sq ft) 


(°F) 

(OF/ft 1 / 2 ) 
(OF/ft 1 / 2 ) .. 
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Imaginary unit, I 


.,1/3 _ 

1 J -l/3 


r 1/3 j 




1/3 

2/3 


> Bessel functions 


IHq 

%' 

k 

L 

l 

n 

Q 

*i 

^o. 

r 

T 

t 

V 

w 

X 

y 


thermal conductivity of metal, Btu/(hr)(sq. ft)(°F/ft) 

blade length, (ft) 

net spacing in relaxation method 

coordinate perpendicular to both radial and axial 
directions 

temperature residual at net point 

heat -transf er coefficient between hot gas and metals 
Btu/(hr)(sq. ft)(°F) 

heat-transfer coefficient between metal and cooling 
air, Btu/(hr)(sq ft)(°F) 

radius, (ft) (also coordinate in radial direction) 
temperature, (°F) 
average thickness, (ft) 
tip speed, (ft/sec) 

linear dimension on rim common with blades, (ft) 
distance of blade point from blade tip, (ft) 

L - x, (ft) 
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z coordinate in axial direction 

a (2 4l /kt b ) l/2 

3 (2^/ktj.) 1 / 2 3 4 

7 (^rtq^/lcApot) 1 / 2 

0 T m -T a 


Subscripts : 

a cooling air 

b blade 

e effective 

m metal 

r rim (with A denotes area exposed to gas) 

rot rotor 


ANALYSIS 

One -Dimensional Temperature Distribution 

A radial -temperature distribution was obtained by the method 
developed in reference 2 for a turbine with and without cooling 
fins on the rim and with blades having long, thin trailing sections 
and a high degree of twist. The following simplifying assumptions 
were made; 

1. Constant mean values of the area, perimeter, theimal con- 
ductivity, and heat-transfer coefficients were used for each 
section. 


2. No temperature gradients other than radial were considered. 

3. The effect of radiation was negligible. 

4. Cooling on both sides cf the turbine rim was considered. 
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5. For the turbine with cooling fins on the rim, the over -all 
heat-transfer coefficient was taken as 4.4 times the flat-plate 
value of q <J . The value 4.4 was calculated from the fin dimensions 

by the method of reference 10 for fins that were so designed as to 
withstand the estimated stresses. The flat-plate value of q Q was 

calculated by Reynolds analogy (reference 11, p. 162) from values 
of the friction coefficient appropriate to the relative velocity of 
the rotor and the air. 

Expressions for the temperature in the blades, rim, and rotor 
of a turbine are derived in reference 2. These expressions were 
obtained by setting up heat balances for an element of the blade, 
the rim, and the rotor, and by solving the resulting differential 
equations. This process yielded the following three equations: 

Blade, 

T e - T m = C cosh ot» (x— B) (l) 


Rim, 

T m - T a = DJ 0 (ipr) + EiH 0 (ipr) (2) 

Rotor, 

T m - T a = r X /2 Fi l/3 J. l/3 (2/3 r ir 3 / 2 ) + Gi" 1 / 3 ( 2 /3yir 3 / 2 ) 

(3) 

The metal temperatures of the blade, rim, and rotor were found for 
ranges of the variables from x = 0.2116 foot (corresponds to 
r = 0.3750 ft) to x = 0 feet (corresponds to r = 0.5886 ft), 
from r = 0.2775 to r = 0.3750 foot, and from r = 0 to 
r = 0;2775 foot, respectively. 

The six integration constants in equations (l) to (3) may be 

evaluated as follows: From the boundary conditions, 0 for 

or 

r = 0 and = 0 for x = 0, it can be shown that B = G = 0. 

The other four integration constants may be found by solving the 
equations obtained by equating temperatures and heat flows at the 
Junctions of the various sections. These equations were found to 
be 



(0.2775) 1 /%i 1 / 3 J_ 1 / 3 2/37i(0. 2775) 3 / 2 = DJ 0 ip (0.2775) + EIHq ip (0.2775) 
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t r = 0.1388, (ft) 

Ayot ~ 0*1722, (sq f't) 

= 30 Btu/(hr)(sq ft) (°F) (without fins) 
= 132 Btu/(hr)(sq ft)(°F)(with fins) 
a = 27.89, (ft) -1 
8=6, (ft)"l (without fins) 

8 = 12.56, (ft) -1 (with fins) 

7 = 13.7, (ft) -3 / 2 


The solutions of equations (4) to (7), upon substitution of these 
values, gave the following values for C, D, E, and F: 



Without 

fins 


With fins 

c = 

2.59° F 


C = 

4.41° F 

D = 

384.60° 

F 

D = 

35.92° F 

E = 

-71.53° 

F 

E = 

-13,360.00° F 

F = 

880.70° 

f/ ( ft) 1 / 2 

F = 

117.30° F/ (f t)^/ 2 


The radial -temperature distributions through the turbine were then 
obtained (fig. 2) by the use of these calculated integration con- 
stants and equations (l) to (3). 

Special attention was then given to the radial -temperature 
distribution through the turbine blades. Stress -to -rupture data 
for cast S816 high -temperature alloy, obtained from figure 3 in 
reference 13, were combined with the blade-stress distribution, ' 
obtained by the use of equation (3) in reference 9, and allowable 
blade -temperature curves were obtained for various blade -tip speeds 
and various blade lives. The allowable blade-temperature curves 
for 1000-hour blade life are plotted in figure 3 with the blade- 
temperature distributions for the turbine with and without cooling 
fins on the rim. 

In order to determine the effect on blade life of the addi- 
tion of cooling fins on the rim, an allowable blade-temperature- 
distribution curve tangent to, or Just above, the blade- 
temperature-distribution curve for the turbine with rim fins must 
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"be obtained for a fixed blade tip speed (equal to the maximum 
allowable blade-tip speed for the turbine without rim fins) and for * 
a blade life as yet unknown. A family of allowable blade-temperature- 
distribution curves was found for this fixed tip speed and for dif- 
ferent blade lives. The curves for 1000-, 2000-, and 3000-hour 
blade life are plotted in figure 4 with the blade -temperature dis- 
tributions for the turbine with and without rim fins at a tip speed 
of 1170 feet per second. 

The effect of the addition of cooling fins on the rim on the 
effective gas temperature was then investigated. In order to deter- 
mine this effect, an allowable blade -temperature curve was fixed for 
1000 -hour blade life and a blade tip speed of 1170 feet per second. 

A family of blade-temperature-distribution curves was then found for 
various effective gas temperatures. A blade-temperature -distribution 
curve tangent to, or just below, the fixed allowable blade rtemperature- 
distribution curve was required. The curves so obtained are shown 
in figure 5. 


Two-Dimensional Temperature Distribution 


Two-dimensional temperature distributions were obtained through 
the turbine rim and rotor by considering temperature gradients in 
both the radial and the axial directions. Symmetry was assumed 
about the radius of the rotor because cooling was considered at 
both sides of the rim. From a heat balance for an element 2«rdrdz 
of the wheel, the heat -flow equation was found to be 


d 2 e 1 be 

ar 2 + r 5? + £ Z 2 


0 


( 8 ) 


Along the cooled edge of the turbine, the boundary condition was 


be 



The boundary condition at the rim surface exposed to the hot gas 
was 


5r " k 


And the boundary condition at the Junction of the rim and blades 
was 
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3 °kAb |f = 2nkr r w |f - q^e 


where account is taken of heat entering the rim by both convection 
and conduction. It was assumed there was no heat flow across the 
other boundaries. 

A solution of equation (8) was found by applying the relaxa- 
tion method (reference 14). The one -dimensional temperature dis- 
tribution previously obtained was assumed as the radial -temperature 
distribution along the cooled edge of the turbine. From these 
values, the definition of 0, and the boundary conditions, esti- 
mated values of 9 were obtained and inserted at each point of a 
square network of points. Temperature residuals (interior heat 
sinks) were then calculated at each net point and reduced to mini- 
mum values by the use of fixed relaxation patterns. From these 
revised values of 9, the metal temperature was found at 

each net point, and isotherms were then drawn. (See fig. 6.) 


Three-Dimensional Temperature Distribution 


In order to obtain a more accurate temperature distribution 
around a blade root, a detailed three-dimensional relaxation was 
carried out for a small portion of the turbine rim. The portion 
selected extended from the center of a blade to a point midway 
between that blade and the next blade, from the center of the rim 
to the cooled edge of the rim, and from the rim surface l/4 inch 
in a radial direction. (See fig. 7.) As before, cooling was 
considered on both sides of the turbine rim. Symmetry was assumed 
about the center of the rim, about the center of the blade, and 
about the line midway between the two blades. Because only a 
small portion of the rim was considered and the curvature of the 
various boundary surfaces of such a portion was slight, the curva- 
ture was disregarded. The heat-flow equation, derived by setting 
up a heat balance for the element cix Sy dz was found to be 


a 2 e a 2 e 

ctoc 2 dy 2 



Boundary conditions were the same as those considered in the two- 
dimensional investigation and the solution was obtained in exactly 
the same way, except that a cubical network of points was required 
here; hence, a new relaxation pattern was employed. 
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RESULTS AND DISCUSSION 

Three different sections of, the specially designed turbine 
blade considered in this investigation are shown in figure 1. The 
root section, the central section, and the tip section are presented 
in order to show the high degree of twist and the long, thin trail- 
ing section of the blade, both of which preclude direct cooling. 

The radial -temperature distribution through the turbine rotor, 
rim, and blades, obtained for the turbine with and without cooling 
fins on the rim by application of equations (l) to (3), is shown 
in figure 2. The addition of fins lowered the rotor and the rim 
temperatures by at least 325°. F, but the reduction in blade tem- 
peratures decreased as the blade tip was approached and, for the 
upper half of the blade, the temperature reduction due to the addi- 
tion of cooling fins on the rim was negligible. 

The temperature distribution through the blade for the turbine 
with and without cooling fins on the rim plotted with allowable 
blade -temperature curves for various tip speeds and 1000-hour blade 
life is presented in figure 3. The critical blade point for the 
turbine without fins and for an effective gas temperature of 1500° F 
was found to be at 0.35 the distance from the blade root; at this 
point, the allowable, blade temperature was about 1435° F and the 
tip speed was 1170 feet per second. The critical blade point for 
the turbine with fins occurred slightly farther from the blade root. 

The effect on blade life of the addition of cooling fins on the 
rim is shown in figure 4. For a fixed blade tip speed of 1170 feet 
per second and an affective gas temperature of 1500° F, the addition 
of fins indicates cui Increase in blade life from 1000 to 2000 hours. 

In practice, such an increase would be difficult to realize because 
the effective gas temperature and weight flow would have to be kept 
uniform and the blade material would have to be homogeneous to a high 
degree. The results indicate that near the rupture point the life of 
the material is highly sensitive to small temperature changes. 

The effect of the addition of cooling fins on the rim on the 
effective gas temperature is illustrated in figure 5. It was found 
that for a blade life of 1000 hours and a blade tip speed of 1170 feet 
per second, the addition of fins permitted an effective gas temperature 
of 1530° F, an increase of 30° F. It is pointed out that figure 6 in 
reference 3, which shows possible increases in effective gas tempera- 
ture for a range of blade parameters, gives in this case a result con- 
sistent with this value. 

From the preceding results, it may be observed that the addition of 
cooling fins on the rim resulted in doubling the blade life, although 
it permitted only a slight increase in effective gets temperature. 
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The results of a two-dimensional investigation of the tempera- 
ture distribution through the turbine rim and rotor, obtained by 
use of the relaxation method, are shown in figure 6. The inves- 
tigation was made primarily to determine the accuracy of the radial 
(one -dimensional) distribution previously found. The results are 
limited to the turbine rim and rotor because no significant changes 
could be found in the blade -temperature distribution without a 
considerable amount of calculation necessitated by a reduction in 
net spacing and a large increase in the number of net points used. 
The radial -temperature distribution proved to be accurate; the 
addition of cooling fins on the rim resulted in lowering the tem- 
perature at the blade root by about 300° F. Correspondingly 
larger increases were found nearer the turbine axis. 

Although the three-dimensional temperature investigation 
through a section of the turbine rim near a blade root is not 
directly related to the other investigations reported herein, it 
is included because of its value for any future study of rim 
thermal stresses. The results of the investigation are shown by 
the isothermal surface in figure 7; the l/4-inch section of the 
rim that was considered is divided into three layers in order to 
show more clearly the isothermal surfaces. A careful study of 
figure 7 shows the temperature changes in the various directions 
for the rim section that was used in the investigation. 


SUMMARY GF RESULTS 

The following results were obtained from a theoretical inves- 
tigation of the temperature distribution through a turbine with a 
specially designed blade and with and without cooling fins on the 
rim: 


1. For a fixed tip speed of 1170 feet per second gnri an effec- 
tive gas temperature of 1500° F, the addition of rim fins permitted 
a sllgto increase in the effective gas temperature (from 1500° to 
1530 F ) * 

2 - The two-dimensional temperature distribution, obtained by the 
relaxation method, proved the analytically determined radial distri- 
bution to be sufficiently accurate for most applications. 

3. The three-dimensional temperature distribution revealed the 
conditions existing in a section of the rim near a blade root and may 
be applied to investigations of thermal stresses. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 



12 


NACA EM No. E7I22a 


REFERENCES 

1. Reuter, J. George, and Gazley, Carl, Jr.: Computed Temperature 

Distribution and Cooling of Solid Gas-Turbine Blades. NACA 
RM No . E7Bllh, 1947 . 

2. Brown, W. Byron: Cooling of Gas Turbines. I - Effects of Addi- 

tion of Fins to Blade Tips and Rotor, Admission of Cooling Air 
through Part of Nozzles, and Change in Thermal Conductivity of 
Turbine Components. NACA RM No. E7Blla, 1947. 

3. Wolfenstein, Lincoln, Meyer, Gene L., and McCarthy, John S. : 

Cooling of Gas Turbines. II - Effectiveness of Rim Cooling of 
Blades NACA RM No. E7Bllb, 1947. 

4. Sanders, J. C-, and Mendelson, Alexander; Theoretical Evaluation 

of Methods of Cooling the Blades of Gas Turbines. NACA 
RM No . E7Bllg, 1947 . 

5. Brown, W. Byron, and Livingood, John N. B.: Cooling of Gas 

Turbines. Ill - Analysis of Rotor and Blade Temperatures in 
Liquid-Cooled Gas Turbines. NACA RM No. E7Bllc, 1947. 

6. Brown, W. Byron, and Monroe, William R. : Cooling of Gas Turbines. 

IV - Calculated Temperature Distribution in the Trailing Part 
of a Turbine B]ade Using Direct Liquid Cooling. NACA 
RM No. E7Blld, 1947. 

7. Wolfenstein, Lincoln, Maxwell, Robert L., and McCarthy, John S.; 

Cooling of Gas Turbines. V - Effectiveness of Air Cooling 
of Hollow Blades. NACA RM No. E7Blle, 1947. 

8. Livingood, John N. B., and Sams, Eldon W. : Cooling of Gas 

Turbines. VI - Computed Temperature Distribution through 
Cross Section of Water-Cooled Turbine Blade. NACA 
RM No. E7BlIf , 1947. 

9. Bressman, Joseph R., and Livingood, John N. B. : Cooling of Gas 

Turbines. VII - Effectiveness of Air Cooling of Hollow Turbine 
Blades with Inserts. NACA RM No. E7G30, 1947. 

10. Biennann, Arnold E., and Pinkel, Benjamin: Heat Transfer from 

Finned Metal Cylinders in an Air Stream. NACA Rep. No. 488 
1934. ' 

11. McAdams, William H. : Heat Transmission. McGraw-Hill Book Co., 

Inc., 2d ed., 1942, pp. 162, 168, 236. 



MCA EM Wo. E7I22a 


13 


12. Jahnke, Eugen, and Emde, Fritz: Funktionentaf eln mit Foimeln 

und Kurven. (Tables of Functions with Formulae and Curves). 
Sec. 8, 3d rev. ed., B. G. Teubner (Leipzig-Berlin), G. E. 
Stechert and Co. (New York), 1938, pp. 226, 235, 236. 

13. Machlin, E. S., and Nowick, A. S.: Stress Eupture of Heat- 

Eesisting Alloys as a Bate Process. MCA TN No. 1126, 1946. 

14. Bnmons, Howard W.: The Numerical Solution of Partial Differ- 

ential Equations . Quarterly Appl. Math.) vol II, no. 3, 

Oct. 1944, pp. 173-195. 




16 


NACA RM No. E7 I 22a 



0123456 7 


Radius, r, in. 



Figure 2. - Rad i al -tempe ratu re distribution through turbine with special 
turbine blades for effective gas temperature of 1500° F. 




Figure 3. - Radi al-temperature distribution through special turbine 
blade. for effective gas temperature of 1500° F. 
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Blade position, y/l 



Figure 5. - Allowable increase in effective gas temperature for turbine 
with cooling fins on the rim. 8lade tip speed, 1170 feet per second; 
blade life, 1000 hours. 
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Figure 7. - Isothermal surfaces in section of turbine rim near blade 
root considered for three-dimensional relaxation. (Temperature 
in ° F.) 



